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Abstract 

We discuss two alternative methods, based on the Lindstedt-Poincare technique, for the removal 
of secular terms from the equations of perturbation theory. We calculate the period of an anhar- 
monic oscillator by means of both approaches and show that one of them is more accurate for all 
values of the coupling constant. 
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I. INTRODUCTION 



Straightforward application of perturbation theory to periodic nonlinear motion gives rise 
to secular terms that increase in time in spite of the fact that the trajectory of the motion is 
known to be bounded^. One of the approaches commonly used to remove those unwanted 
secular terms is the method of Lindstedt-Poincare^, recently improved by Amore et al££ 
by means of the delta expansion and the principle of minimal sensitivity. 

There is also another technique that resembles the method of Lindstedt-Poincare which 
is suitable for the removal of secular terms^. Discussion and comparison of such alternative 
approaches may be most fruitful for teaching perturbation theory in advanced undergraduate 
courses on classical mechanics. 

In Sec. |n]we present a simple nonlinear model to which we apply the alternative pertur- 
bation approaches in subsequent sections. In Sec. IHII we apply straightforward perturbation 
theory and illustrate the outcome of secular terms. In Sec. IIVI we show how to remove 
those secular terms by means of the Lindstedt-Poincare method 1 '^. In Sec. we develop 
an alternative method that appears in another textbook^ and that closely resembles the 
method of Lindstedt-Poincare. In Sec. IV II we describe an improvement to the method of 
Lindstedt-Poincare proposed by Amore et al££. Finally, in Sec. IVIII we compare the period 
of the motion calculated by all those approaches. 

II. THE MODEL 

In order to discuss and compare the alternative perturbation approaches mentioned above, 
we consider the simple, nonlinear equation of motion 

x(t) + x{t) = -Xx 3 (t) (1) 

with the initial conditions x(0) = 1 and x(0) = 0. In Appendix 1X1 we show that we can 
derive this differential equation from the equation of motion for a particle of mass m in a 
polynomial anharmonic potential with arbitrary quadratic and quartic terms. Notice that 
E = x 2 /2 + x 2 /2 + Ax 4 /4 = 1/2 + A/4 is an integral of the motion for (JTJ) and that the 
motion is periodic for all A > —1 for the initial conditions indicated above. 
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III. SECULAR PERTURBATION THEORY 



The straightforward expansion of x(t) in powers of A 



x(t) = -£ Xj (t)\ j (2) 

3=0 

leads to the perturbation equations 

x (t) + x (t) = 
xi(t) +x 1 (t) = -xl(t) 

n— 1 j 

x n (t)+x n (t) = - ^2 ^2x k (t)xj- k {t)x n -j-i(t), n = 2,3, . . . (3) 

j=0 k=0 

with the boundary conditions Xj(0) = 5jq and Xj(0) = for all j > 0. Clearly, the solution 
of order zero is x (t) = cos(i). 

All those perturbation equations are of the form y(t) + y(t) = f(t), where fit) is a 
linear combination of cos(jt), j = 1,2,.... Such differential equations, which are commonly 
discussed in introductory calculus courses, are in fact suitable for illustrating the advantage 
of using available computer algebra systems. 

It is well known that any term proportional to cos(t) in f(t) gives rise to a secular term 
after integration^. For example, since x^(t) = [cos(3t) + 3cos(t)]/4 we obtain 



/ x 1 
Xl{t) = - 



cos(3t) — cos(t) . , 

— - 3t sin(t) 



(4) 
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that clearly shows that grows unboundedly with time in spite of the fact that the 

exact motion is periodic for all A > —1. 

IV. METHOD OF LINDSTEDT POINCARE 

There are several suitable mathematical techniques that overcome the problem of secular 
terms mentioned above^. For example, the method of Lindstedt-Poincare is based on the 
change of the time variable 

s = (5) 

where ^7 plays the role of the frequency of the motion, and, therefore, the period T is given 
by 

T = — ■ (6) 

x/7 



The equation of motion (JTJ thus becomes 

7x"(s) +x{s) = -Xx\s) (7) 

where the prime stands for differentiation with respect to s. If we expand both x(s) and 7 
in powers of A 

00 

7 = E7 3 A J (8) 

j=0 

with 70 = 1, then we obtain the set of equations 
x'^s)+x (s) = 

= -zjj(s) -71X0(5) 

n— 1 J n 

<(«)+a ; n(s) = -S£^fc(s)»i-fc(s)a;n-j--i(s)-Z]77^-j(s), 

j=0 fc=0 j=l 

n = 2,3,... (9) 

We choose the value of the coefficient 7„ in order to remove the secular term from the per- 
turbation equation of order n. For example, it follows from — Xg — 7iXq = (71 — 3/4) cos(t) — 
cos(3t/4) that 71 = 3/4 is the right choice at first order. Proceeding exactly in the same 
way at higher orders we obtain the coefficients 

71 = 4 , 7 2 = -j^ , 73 = ^ (10) 



and the approximate period 



- , 32 ^ (id 

; (9A 3 - 12A 2 + 384A + 512) 



V. ALTERNATIVE LINDSTEDT POINCARE TECHNIQUE 

Perturbation theory provides a A-power series for the frequency of the motion for 
example, for our model it reads 

u 2 = l + w 1 X + w 2 X 2 + ... (12) 

An alternative perturbation approach free from secular terms is based on the substitution 
of this expansion into the equation of motion $1} followed by an expansion of the resulting 
equation 

x(t) + (u 2 - wi\ - w 2 X 2 + ...) x(t) = -Xx 3 (t) (13) 



in powers of A, as if u were independent of the perturbation parameter 5 . The perturbation 
equations thus produced read 

x (t) + co 2 x (t) = 

xi(t) + uj 2 xi(t) = —xl(t) + w 1 x (t) 

n—1 j n 
X n (t) + U 2 X n (t) = - XkfyXj-kfyXn-j-!^) - ^2wjX n -j(t), 

j=0 k=0 j=l 

n=2,3,.... (14) 

Notice that xq = cos(ut) depends on u and so does each coefficient w n that we set to remove 
the secular term at order n. Consequently, we have to solve the partial sums arising from 
truncation of the series (fT^j) for u in order to obtain the frequency and the period 

T=~. (15) 

in terms of A^. 

A straightforward calculation through third order yields 

m = I , «* = -i2jbi ' W3 = (16) 



from which we obtain 



f V30A 2 + 96A + 64 + 2(3A + 4) 
- = ^ - A (IT) 



VI. VARIATIONAL LINDSTEDT POINCARE 

Amore et al£& have recently proposed a variational method for improving the Lindstedt- 
Poincare technique. It consists of rewriting equation (JTJ) as 

x{t) + (l + « 2 ) x(t) = 5 (-Ax 3 (t) + a 2 x(t)) (18) 

where a is an adjustable variational parameter, and 5 is a dummy perturbation parameter 
that we set equal to unity at the end of the calculation. When 6 = 1 the modified equation 
of motion (|T8|) reduces to equation (JTJ that is independent of a. Following the Lindstedt- 
Poincare technique we change the time variable according to equation (|SJ) thus obtaining 

7x"(s) + (1 + a 2 )x(s) = 5 (-Ax 3 (s) + a 2 x(s)) . (19) 



We then expand both x and 7 in powers of S and proceed exactly as is Sec. IIV| except that 
in this case 70 = 1 + a 2 . Thus we obtain 

4(s) + x (s) = 

x"(s) + = — (-Xxq(s) - Jix'q(s) + a 2 x (s)) 

7o v ' 

1 / n_1 J 

x" n {s) + x n (s) = — a;fc(s)x j _ fc (s)x„_ j _i(s) 

j=0 k=0 



To 



-EKj(«) + «Vi(s) , 7i = 2,3,.... (20) 
i=i / 

The period of the motion is given by equation (JHJ). 

Choosing the value of 7 n in order to remove the secular term from the perturbation 

equation of order n we obtain 

3A - 4a 2 3A 2 3A 2 (3A - 4a 2 ) 

71 ~ I ' 72 ~ 128 (l + « 2 ) ' 73 " 512 (1 + a 2 ) 2 ' { ' 

Notice that these coefficients reduce to those in equation (jlOj) (multiplied by the proper 

power of A) when a = 0. Since the actual value of 7 is independent of a when 5 = 1, we 

make use of the principle of minimal sensitivity developed in Appendix [BJ The root of 

^X>, = (22) 

is 



and we thus obtain 



For this particular problem we find that the value of a given by the PMS condition 
is such that 72^+1 = for all j > 0. 



« = ^ (23) 

T = 8^tt(3A + 4) 

/ (207A 3 + 852A 2 + 1152A + 512) ' 



VII. RESULTS AND DISCUSSION 



Fig. ^ shows the period as a function of A given by the perturbation approaches discussed 
above and by the exact expression 1 



We appreciate that the variational Lindstedt-Poincare method^ yields more accurate re- 
sults than the straightforward Lindstedt-Poincare technique^, and that the alternative 
Lindstedt-Poincare approach 5 - is the best approach, at least at third order of perturbation 
theory. 

All those expressions yield the correct value T(0) = 2tt and become less accurate as A 
increases. However, two of them give reasonable results even in the limit A — > oo. The exact 
value is 

lim ^\T exact = 2V2 / . = 7.4162987. (26) 

A ^°° Jo ./1 + CO s(0) 2 



cos(fl) 

The standard Lindstedt-Poincare technique fails completely as shown by 



lim V\T LP = 0. (27) 



A^oo 



The variational improvement proposed by Amore et al£& corrects this anomalous behavior 



lim y/XT VLP = = 7.4112410. (28) 

A^oo 23 

Finally, the alternative Lindstedt-Poincare method 5 gives the closest approach 



lim V\T ALP = Ti 

A^oo 



\ 



64 _32V30\ =741g5905 (2Q) 



APPENDIX A: DIMENSIONLESS EQUATIONS 

Transforming an equation of physics into a dimensionless mathematical equation is most 
convenient for at least two reasons. First, the latter is much simpler and reveals more clearly 
how it can be solved. Second, the dimensionless equation exhibits the actual dependence of 
the solution on the parameters of the physical model. 

In order to illustrate how to convert a given equation into a dimensionless one we consider 
a particle of mass m moving in the potential 

v(q) = f g 2 + ^s 4 . (Al) 

The equation of motion is 

mq = —V2q — v^q 3 (A2) 
and we assume that q(0) = q and q(0) = v . 



We define a new independent variable s = uot + <fi, where is a phase, and ojq = \jvijm 
is the frequency of the motion when t> 4 = 0. Suppose that q = and q — A at t — tx, 
then we define the dependent variable x(s) = q(t)/A and choose = — c^o^i so that x(s) is 
a solution of the differential equation 

x"{s) +x{s) = -Xx{s) 3 (A3) 

where A = V4,A 2 /v2 and the initial conditions become x(0) = 1 and x'(0) = 0. 

The dimensionless differential equation (|A3J) resembles the equation of motion for a par- 
ticle of unit mass moving in the potential V(x) = x 2 /2 + Ax 4 /4. Its period T depends on 
A and, therefore, the expression for the period of the original problem T' = T/ujq clearly 
reveals the way it depends upon the model parameters m, v%, v^, and A. 



APPENDIX B: VARIATIONAL PERTURBATION THEORY 

Variational perturbation theory is a well-known technique for obtaining an approximation 
to a property -P(A) in a wide range of values of the parameter A. Suppose that -P(A) is a 
solution of a given equation of physics F(X, P) = that we are unable to solve exactly. In 
some cases we can obtain an approximation to -P(A) in the form of a power series -P(A) = 
P0 + -P1A + . . . by means of perturbation theory. If this series is divergent or slowly convergent 
we may try and improve the results by means of a resummation technique. 

Variational perturbation theory consists of modifying the physical equation in the form 
_F(£, a, A, P) = 0, where a is a variational parameter (or a set of them in a more general 
case) and £ is a dummy perturbation parameter so that F(l, a, A, P) = F(X, P). 

Then we apply perturbation theory in the usual way, calculate N + 1 coefficients of the 
perturbation series 

00 

P(£,e*,A) = 5»,A)e (Bl) 
3=0 

and construct an approximation of order N to the property 

N 

pW(a,\) = E P j(®^)- (B2) 

j=Q 

If the partial sums P^(a,X) converged toward the actual property as N — > 00, then 
pt°°'(a, A) = -P(A) would be independent of a. However, for finite N the partial sums 



S 



do depend on the variational parameter a. It is therefore reasonable to assume that the 
optimum value of this parameter should be given by the principle of minimal sensitivity 
(PMS)£: 

ft 

= 0. (B3) 



Ap M (a,A) 



a=a N (\) 

In many cases P^^a^A), A) converges towards -P(A) as Af -> oo, and, besides, 
P'^(aiv(A), A) behaves like P(X) with respect to A even at relatively small perturbation 
orders. 
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FIG. 1: Absolute error on the period \T exac t — T approx \ as a function of the coupling A for the 
Lindstedt-Poincare (LP, solid thin line), the Alternative Lindstedt-Poincare (ALP, solid bold line) 
and the Variational Lindstedt-Poincare (VLP, dashed line) methods. 
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